Spin-triplet s-wave local pairing induced by Hund's rule coupling 



J. E. Han 

Department of Physics, The Pennsylvania State University, University Park, PA 16802-6300, USA 
Department of Physics, State University of New York at Buffalo, Buffalo, NY 14260, USA 

(Dated: February 2, 2008) 

We show within the dynamical mean field theory that local multiplet interactions such as Hund's 
rule coupling produce local pairing superconductivity in the strongly correlated regime. Spin-triplet 
superconductivity driven by the Hund's rule coupling emerges from the pairing mediated by local 
fluctuations in pair exchange. In contrast to the conventional spin-triplet theories, the local orbital 
degrees of freedom has the anti-symmetric part of the exchange symmetry, leaving the spatial part 
as fully gapped and symmetric s-wave. 
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I. INTRODUCTION 

For the last two decades, discovery of unconventional 
superconductors have fueled intensive research for novel 
pairing mechanisms. Of particular interest have been 
the cupratei and heavy-fermion^ systems where the very 
existence of superconductivity is quite surprising due to 
the strong Coulomb repulsion within the framework of 
the conventional Eliashberg-Migdal theorji^. Recently 
discovered class of superconductors close to the ferro- 
magnetic phase such as Sr2RuO^, UGe^^ demand better 
understanding between symmetry and superconductivity 
which goes beyond the s-wave BCS theory. 

One of the most puzzling aspects of the supercon- 
ductivity in the strongly correlated regime is the role 
of the repulsive Coulomb interaction. Common belief 
is that the pairing mechanism is driven or assisted by 
the Coulomb repulsion to an attractive interaction be- 
tween electrons which results in spatial or temporal struc- 
tures in the superconducting wavefunction to avoid di- 
rect Coulomb repulsion. Even some conventional s-wave 
phonon-mediated superconductors such as alkali-doped 
fuUerenes have also shown extraordinarily high transition 
temperatures^ despite the strong Coulomb repulsion near 
the Mott transitioniiS. 

Even before the onset of superconductivity, direct 
Coulomb interaction is heavily renormalized away by 
electrons moving in a correlated manner to form a narrow 
resonance near chemical potential. The superconductiv- 
ity in such regime should be understood in terms of those 
correlated electronic basis, typically called quasi-particles 
(QP). The problem of mutual renormalization effects be- 
tween electron-electron and electron-phonon (el-ph) has 
received a great deal of interest resulting in different in- 
terpretations, ranging from the conventional Coulomb 
pseudopotential^, field-theoretic approachesiSiii to nu- 
merical analysisi^. Recently, it has been pointed oulii^ 
that the symmetry of the el-ph interaction critically influ- 
ences the pairing interaction where the two interactions 
are so interwoven that the effective pairing interaction 
cannot be reduced into a form similar to the McMillan 
formulaM. 

We show here on the basis of generic hamiltonians that 



local orbital correlations are crucial for the superconduc- 
tivity and, in particular, predict the Hund's rule (HR) 
induced superconductivity. We discuss the local pairing 
as mediated by self-generated fluctuations within the on- 
site interactions and show that the Coulomb interaction 
and multiplet interaction work in fundamentally differ- 
ent way than described in perturbative theories. The HR 
superconductivity has the spin-triplet pairing state. In 
contrast to the well-known non-s-wave triplet pairing^-^, 
the anti-symmetric part of the exchange symmetry of the 
pairing wavefunction is taken up by the local orbital vari- 
ables, leaving the spatial part as the symmetric s-wave. 
We will discuss possible consequences of this new pairing 
structure. 

This paper is organized as follows. In sectional we 
define the problem in models of multiplet interactions 
for Hund's rule and Jahn- Teller el-ph couplings. In the 
following section ITTll the pairing mechanism is explained 
in a perturbation theory of a simplified dynamical mean 
field theory (DMFT). Results from full DMFT calcula- 
tions are presented in the section llVl and are discussed in 
detail. In Appendices IXknd IbI we describe the calcula- 
tions for the local energy levels and perturbation theory. 
Finally, a new Hubbard-Stratonovich decoupling scheme 
for the Hund's rule coupling which bypasses the sign- 
problem is introduced in Appendix ICl 

II. MODELS 

We consider cases of doubly degenerate band electrons 
at half-filling which are locally coupled via Hund's rule 
(HR) and Jahn- Teller (JT) phonons. We model the 
multi-band electronic system with the on-site Coulomb 
repulsion of strength U as 

ijma i 
/ i 

where dj^^ (dima) is the electron creation (annihilation) 
operator acting on site i, orbital m (= 1, 2), spin a, fj, is 



the chemical potential and Ni = X^mo- '^j'"<^- "^^^ ^^P" 
ping integral tij is chosen to give a semi-elliptical density 
of states with the half-bandwidth W/2 = The local 
multiplet interaction is taken as Hmuit ^^hr or i/j-] 
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where the site index i has been omitted for brevity. K 
is the HR coupling constant and Lp^ [v = 1,2) are the 
phonon fields with the bare frequency ojph, and g the 
el-ph coupling constant. In the anti-adiabatic limit of 
phonons (wph — > oo)^^, iTjT maps to the same form 
as -f^HR but with a negative (fictitious) Kjt = — 2A 
(A = f/^/wph). Therefore, the el-ph superconductivity 
is suppressed as the physical HR coupling is turned up. 
We will show that, as K increases further, superconduc- 
tivity re-emerges in a spin-triplet channel. 



III. LOCAL PAIRING MECHANISM 

We describe main ideas behind the local pairing before 
presenting results from the full dynamical mean field the- 
ory. In the strongly correlated regime, QP states within 
the energy band [—zW,zW] {z the QP renormalization 
factor, W non-interacting bandwidth) are usually respon- 
sible for low-energy manybody phenomena. The effect of 
U is implicit in the renormalization factor z and the trun- 
cated high energy states. In such limit, we simplify the 
DMFT by singling out the QP band and ignore incoher- 
ent high energy excitations as 
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where cj^^^^ (Cmko-) is creation (annihilation) of a QP 
state, is QP energy (absorbing chemical potential), 
tq zW) is hopping integral on and off the impurity 
and Ns is the number of sites. Contributions from the 
high energy states are correctly considered in the full 
DMFT calculations. Hi^t is the interacting part of the 
original hamiltonian. Note that the interaction only acts 
on the impurity site. For small QP bandwidths with zW 
smaller than other energy scales, we further simplify the 
above effective hamiltonian to a two-site model. 
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FIG. 1: Pair correlation energy Ep for two-site model 
Eq. with the Hund's rule coupling. Ep is attractive for 
a wide range of hopping integral to the quasi-particle bath, 
tq. Coupling to quasi-particle bath (thick line) produced a 
good agreement with the perturbation theory, —30tq/KU*^ 
{U* = U + 2K). The optimal interaction is at 4*^/(7 ~ K for 
the maximum pair attraction Ep oc —K. 
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the QP energy (e, « 0), where the narrow QP band is 
represented by a single energy level. H^g has the simi- 
lar construction as the exact diagonalization implemen- 
tation of the DMFT, where the bath spectrum is sampled 
with a finite number of energy levelsii. Representing the 
QP states with several discrete levels over a small energy 
range made little change. This procedure is expected to 
be valid as long as the electrons slow down by the strong 
Coulomb interaction so that the local picture is meaning- 
ful and there is a QP reservoir which supplies particles 
at an infinitesimal energy cost. 

We solve the two-site problem Eq. Q using the exact 
diagonalization method. We analyze the solution as a 
perturbation expansion in tc^. We estimate the pairing 
correlation by computing the energy gain Ep for creating 
two extra electrons on a half-filled impurity site. 



p = [E{Nd + 2) - EiNd)] - 2[E{Nd+l) - E{Nd)], (6) 



where E{n) is the ground state energy with n total elec- 
trons and Nd is the orbital degeneracy {Nd = 2). It 
is crucial to have multiple electrons per site for the lo- 
cal multiplets to induce the local pairing. As shown in 
Refii^, the local pairing has a strong filling dependency. 

Numerical results shown for the HR problem are with 
parameters of i7 = 4 and K = 0.1. For small tq in a HR 
system, Ep is negative (thick line in Fig.^l, i.e., it is en- 
ergetically favorable to create mutually interacting quasi- 
particles than to create two independent quasi-particles. 
The impurity ground states are always predominantly of 
filling Nd{= 2) with the total spin S = 1 since any ex- 
tra electrons created on the impurity are repelled by the 
Coulomb interaction and are immediately dissolved into 
the QP bath. Therefore we first examine the atomic limit 
{tq = 0) . Without i/muit the degeneracy of two-electron 
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FIG. 2: Single site energy level scheme for two on-site elec- 
trons {Nd = 2) with Hund's rule and Jahn- Teller couplings. 
The degeneracy of 6 is lifted by the multiplet interaction 
^^muit- The ground states are spin-triplet (5* = 1) and spin- 
singlet {S = 0) for Hund's rule. The Jahn- Teller multiplets 
have inverted level structure of the Hund's rule coupling. 




S=ff~ 
intermediate state 

FIG. 3: One of the contributions in E{Nd + 2) for the pair cor- 
relation energy Ep in Eq. 12UII . Two extra electrons created 
on the interacting site dissolve into the quasi-particle bath to 
avoid Coulomb repulsion, which produces the ground state 
configurations shown on the left and right panel. The two 
unperturbed degenerate ground states couple via exchange of 
the spin-triplet electron pairs in the 4-th order perturbation 
of tq (from the left configuration to the right). In the inter- 
mediate state, the spins on the interacting site are in a locally 
excited state (See Fig. I^J. 



configurations is 6. The HR coupling lifts the degener- 
acy with the ground states of spin-triplet {S = 1) on the 
impurity site with 



-t-rn) — ^ [i]^|U2T' V IT 2i "2T li J ' "li 



Ivac) 



(7) 

with m = +1,0,-1, respectively (see FIG. |21 and Ap- 
pendix 0. The 1st excitation energy is 2K for the 
HR coupling. The local energetics is described in Ap- 
pendix ^ Upon adding two electrons to the interacting 
site, the extra charges are immediately pushed to the QP 
reservoir to avoid the Coulomb repulsion. Through inter- 
site hopping, the ground state triplet degeneracy is lifted. 
The resulting leading order ground state wavefunction 
has the most symmetric form {U ^ K ^ tq > 0), 
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\To) ^ \to) + \Ti) <g> , 
(8) 

where the second ket-vectors for the QP states in the 
direct products have similar definitions as Eq. (|7|). We 
analyze the perturbed ground state wavefunction and its 
energy in terms of the perturbation theory developed in 
Appendix The first order perturbed wavefunction be- 
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with the ionization energy U* /2 = U/2 + K. This in- 
termediate state has an ionized impurity at the excited 
energy U* . With Eq. (|BTT1) and Eq - Hq = -U* /2, one 
gets the second order contribution to the energy 
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where the figure 8 is from the number of terms in Eq. ||^. 
The second order correction to the wavefunction from 
Eq. I|B10|I becomes 
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where \Sa), \Sb) and \Sc) are spin-singlet excited states 
on the interacting site as defined in Appendix IXI The QP 
configurations \sa), \sb), |sc) are defined with the same 
symmetry as the impurity states. 

One of the processes in the second order contribution is 
depicted in Fig. 13 Typical (unperturbed) basis states in 
the ground state are spin-triplets as shown on the left and 
right side. The two states couple via the electron hop- 
ping where a second-order intermediate state is shown in 
the middle. After two hoppings {eg. the processes 1 and 
2 from the left configuration), the resulting intermedi- 
ate state becomes a charge- neutral spin-singlet {S — 0), 
which correspond to the term |S'a) ® |sa) in Eq. Hll|l and 
a local transition between the multiplets in Fig. |21 Other 
charge-excited intermediate states have small contribu- 
tions to the transition amplitude (of order K/U smaller). 

From Eq. ljB12p . the forth order correction to the en- 
ergy becomes 
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for K -^U* . The ground state energy E{Nd+2) becomes 



E{Nd + 2) = U- 
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The first term is the (unperturbed) atomic- limit. The 
last term contributes to the pairing interaction. The 
denominator in Eq. (|13|l comes from two ionized states 
(with charging energy U* after the hoppings 1 and 3 in 
Fig. 131 for example) and one HR-excited state (with exci- 
tation energy ^ 2K after the hopping 2). 
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Ground state for one extra (spin-up) electron can be 
obtained similarly. For hopping integrals smaller than 
other energy scales, the interacting site is predominantly 
of filling 2 and of spin-triplet. Therefore the ground 
state consists of \To) (Ei Icmt) and \Ti) (g) \cm'i) (|c„ia) = 
cl^^\vac)) . Due to the hopping interaction, states of total 
spins Sz = ±3/2 do not appear in the ground state. An 
explicit calculation gives the unperturbed ground state. 
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U ^ A tq, which retains the similar form as Eq. H20() 
for the HR coupling. Leading order wavefunctions and 
their perturbed energies are written down in Appendix IbI 
On the other hand, for a pure Hubbard model without 
the muhiplet terms (iJmuit ^0),Ep^ MUyU^ > sug- 
gesting no superconductivity, which is also supported by 
full DMFT calculationsi^. With iJ,„uit = 0, the electrons 
only interact with charge fluctuations making the local 
internal fluctuations irrelevant. The above perturbation 
results are summarized in Table 1. 



The larger coefficient for |Ti) (g) |c2x) reflects more paths 
for hopping transitions for this configuration. Similar 
calculations as those with Nd + 2 electrons give 
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The state without extra electrons results in no energy 
gain in the 4-th order. 
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Finally, the pairing correlation Eq. © becomes 
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which is in excellent agreement with the exact diagonal- 
ization results in the tq <^ K limit. It is interesting to 
note that Ep becomes smaller for stronger K due to the 
higher multiplet excitation energy. The exact diagonal- 
ization results show maximum attractive pair correlation 
for 4:tq/U ^ K, where the effective hopping exchange 
energy is comparable to the internal excitation energy. 
Upon substitution for tq in Eq. (I^OJ, this relation gives 
Ep oc —K at the optimal pairing. 

One can also interpret the pairing as arising from pair 
exchange. In Fig. O the spin-up {Sz — 1) pair on the 
interacting site hops out via processes 1 and 3, while 
the spin-down pair hops back in via 2 and 4. The charge 
neutral intermediate state forces the temporal overlap be- 
tween the hoppings of the pairs, producing the attractive 
interaction when viewed as a second order perturbation 
in terms of the pair-basis. The same general reasoning 
should apply to other types of local interactions in the 
strong-t/ limit. 

One can carry out the same perturbation analysis for 
the Jahn-Teller interaction Eq. ||3J| , where the pairing en- 
ergy gain Ep is always negative. Ep = — Tt^/AC/*^ for 



TABLE I: The pair correlation energy for an effective 2-site 
2-band model with Hund's rule (K > 0) and the Jahn-Teller 
el-ph (A > 0) coupling. For large U, both of the couplings 
(l^oduce attractive pairing energy. (A = g^/ii^ph.) 

^)As the pairing correlation energy Eq. (|2U|I and Table 
r suggest, the Coulomb interaction and the multiplet in- 
^j^^action (including the JT el-ph interaction) cannot be 
separated into a form similar to A — /i* as in the McMillan 
formula^^, where the Coulomb interaction represented as 
the Coulomb pseudopotential /i* interact additively with 
the multiplet interaction {eg. el-ph coupling) represented 
by the mass-renormalization factor A. 



IV. PAIRING INSTABILITY: DMFT RESULTS 

Although the 2-site model is suggestive of the nature 
of the pairing, we show the existence of superconduc- 
tivity from a lattice calculation. The lattice hamilto- 
nian Eq. |^ is solved within the DMFT which approx- 
imates one-particle self-energies and vertex corrections 
as momentum-independent and maps a lattice model to 
an effective quantum impurity problem. We solve the 
impurity problem using quantum Monte Carlo (QMC) 
technique without making any physical approximations 
other than discretizing the imaginary time. It has been 
known that the QMC technique often suffers from the 
sign-problemiS. DMFT-impurity models with the Hub- 
bard and electron-phonon couplings have been known not 
to cause any serious sign problem for most of the physi- 
cal parameter spacc^'^. However, inclusion of the Hund's 
rule coupling terms in Eq. Q produced severe sign- 
problems^*^ when each of the interaction terms are decou- 
pled by the discrete Hubbard- Stratonovich transforma- 
tions!, despite that the decoupled HR terms resemble the 
Jahn-Teller coupling. To overcome such problem, a new 
scheme of mixing discrete^^ and continuous^^ Hubbard- 
Stratonovich transformations has been used. This re- 
moved the problem with the average signs larger than 0.9 
for most of the runs. Details of the procedure is given 
in Appendix |21 It is speculated that the simulation is 
more stabilized by the built-in symmetry of orbitals in 
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the new transformation in contrast to the poorly pre- 
served orbital symmetry in incomplete QMC samplings 
of the old scheme. 

The superconducting instability is probed by comput- 
ing uniform pair susceptibility x- X foi' the spin-triplet 
channel can be expressed in a Bethe-Salpeter equation as 



W/, k", In 



X = xo + xorxo + xorxorxo 
= Vxo(/- Vxorvxo)"^ 



(21) 



where we have used the symmetrized formS^. xo is the 
propagator of two electrons independently propagating 
with zero net momentum, and T is the effective pair- 
ing interaction. (See Fig. ^) Eq. (|21|l is a shorthand 
for matrices with indices over the Matsubara frequency 
ujn- For the static uniform susceptibility, we only need 
in-coming electrons of net zero momenta and net zero 
frequency, as in Fig. ^ Due to the momentum and fre- 
quency conservation at each internal vertices of Feyn- 
man diagrams, any internal pair of legs in the ladder 
diagrams in Fig. ^ also have net zero momenta and fre- 
quencies. The uniform (uncorrelated) 2-particle prop- 
agator xo can be computed with the 1-particle Green 
function as xo(w„) = Yl,]iG{\<i,LUn)G{—]i.,—LOn)- The k- 
summation is performed within xo due to the locality of 
the DMFT. F is approximated by a local quantity r'°'^ 
obtained from a local relation similar to Eq. H21|l . 
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Xo"^ is the local (uncorrelated) 2-particle propagator, 
X^S^iuJn) = G'°'=(w„)G''°'=(-w„), and the local triplet pair 
susceptibility x'"*^ is directly computed from QMC as 

X'°=(ri, r2, rg, r,) = {T[d\^{n)dt{T2)d2.{T:,)d,„{Ti)]) . 

(23) 

Corresponding definition of x'°^ for the singlet channel 
can be found in Ref>i^. Since x'°'^ are directly obtained 
from the QMC measurements and xo, Xo"^ are easily 
computed from the 1-particle self-energy, we derive the 
effective interaction kernel F (= r'°'^) and the uniform 
pair susceptibility x from Eqs. (|21I22|) . 

The transition temperatures are determined by when 
the maximum eigenvalue of the superconducting kernel 
\/Xo in Eq. (|?T|l approaches 1. For such condition, 

the effective interaction must be attractive (F > 0). The 
other requirement is the finite weight of quasi-particles 
at the chemical potential as in Xo('^) ~ 'kzN{0)/uj, with 
7V(0) non- interacting density of states at the chemical po- 
tential. As the system becomes strongly correlated, xo 
becomes smaller as z ^ while the effective local corre- 
lation tends to get larger as electrons are more localized. 
The resulting Tc is determined by the balance between 
Xo and F. 

The attractive interaction V^ff(ijj„, Wm) (= 
-(l/T)F(w„,tj„0) is shown in Fig. for U/W = 0.6 
and K/W — 0.1. It is quite interesting that 14ff(w„, Wm) 
resembles that for the electron-phonon coupling in the 
Eliashberg theoryi where T4ff(tj„, Wm) — g^D{ujn — ^m) 




FIG. 4: Expansion of spin-triplet pairing susceptibility x- The 
local index ma (m = 1,2) is for orbital and spin labels. 




FIG. 5: The effective interaction Vcff (t^;™, i^m) (= 
~r{ujri, LL!m)/T). The attractive interaction is for small 
energy transfer {uj„ ~ ujm)- The interaction at higher 
frequencies is repulsive. 



with the phonon propagator D(w„). The attractive part 
is formed in a uniform attractive valley along the small 
energy transfer region a;„ ~ lo„i. For pure Hubbard 
interactions, the effective interaction Vcsi^m^m) shows 
a very different profile with high energy structures24 
away from the low energy scattering a;„ = ujm- The 
attractive interaction along the Un = line clearly 
demonstrates that the HR-induced local pairing is medi- 
ated by the effective low-energy coupling medium which 
has been self-generated by the multiplet interaction. 
Since the local spin excitation energy is 2K (see Fig. ^ 
and the electron-spin coupling constant is of order 
K, a naive estimate for the effective interaction gives 

K^{2/2K) - K{= 0.2), as in the electron-phonon 
theory. The numerically obtained VcS at 2 — 4 is about 
an order of magnitude larger than the simple estimate. 
This suggests that the renormalization effect is very 
strong and is consistent with the above observation that 
the suppressed hopping due to the Coulomb interaction 
enhances the local interaction^^. A thorough exam- 
ination of the renormalization effects is necessary to 
fully understand these numerical results. Diagrammatic 
methods will prove particularly useful when backed up 
by the QMC and the exact diagonalization methods 
since the basic mechanism of the local pairing is given 
as in Section III. 

Fig- Eta) shows transition temperatures Tc for the HR 
coupling as a function of the Coulomb repulsion U where 
the superconductivity becomes stronger with U. The 
HR superconductivity emerged only in the spin-triplet 
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FIG. 6: (a) Superconducting transition temperature 
Coulomb repulsion U with Hund's rule coupling. F( 
Hund's rule (HR) coupled system, Tc monotonically inc± 
with U showing local pairing mechanism. Gray symbol is 
an extrapolated estimate from the superconducting kernel 
XoT. (b) Quasi-particle renormalization factor z (left axis) 
and effective pairing interaction VcS at lowest Matsubara fre- 
quency (right axis), z decreases as U approaches Uc- The 
product zVcS for the HR system increased as ?7 Uc, consis- 
tent with the increasing Tc. W is the non-interacting band- 
width. 
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FIG. 7: (a) Superconducting transition temperature Tc vs 
Coulomb repulsion U with Jahn-Teller electron-phonon cou- 
pling. The Jahn-Teller (JT) phonon coupling produces su- 
perconductivity with a cross-over from the conventional su- 
perconductivity (small U) to the local pairing regime (large 
U) for 2A/W — 0.1. For a larger electron-phonon coupling 
{2A/W = 0.15), the Mott insulator transition happened at 
smaller U. (b) Renormalization factor z as in Fig.|S] 



channel. Tc monotonically increased as U approached 
the critical value Uc/W « 0.75 for K/W = 0.15, a 
strong indication of the local pairing. Tc/W less than 
0.002 (gray symbol) has been estimated by extrapolat- 
ing eigenvalues of the superconducting kernel ^/xo^\/Xo- 
The HR-coupled superconductivity is strong only in the 
local pairing regime, since at small U there exists no ex- 
ternal low-energy pairing medium in the electronic HR 
coupling. 

14fi(7rT, ttT) at the lowest Matsubara frequency at 
T/W = 1/120 (triangles in Fig.E^b), to the right-hand- 
side scale) show a growing local pairing attraction as 
z — > OSSi. The product zVcff (ttT, ttT) (ex xo^) increases 
a,s U ^ Uc, consistent with the increasing Tc- We em- 
phasize that the local pairing becomes apparent even at 
moderately correlated regime of 2; ~ 0.5, although the 
heuristic argument in section UTTl was given for z — > for 
the sake of clarity. 

We can make some predictions based on general prop- 
erties of the local spin-triplet pairing. (1) The exchange 
anti-symmetry is in the orbital space, in contrast to the 
p-wave spin-triplet theorjii^. Due to the discreteness of 
orbitals, the quasi-particle excitation is gapped. The spa- 
tial part has the s-wave symmetry. (2) For large orbital 
degeneracy (Nj^ > 2), coupling between different pairs 
of orbital indices will likely result in multiple gaps. (3) 
The local pairing is strong for high electron/hole dop- 



inggi^. For small dopings, ferromagnetic fluctuation will 
dominate the local triplet pairing. (4) As can be inferred 
from Fig|21 strong ferromagnetic order suppress the pair- 
exchange. Therefore, if the pairing and ferromagnetic 
phases ever co-exist, the superconducting phase should 
be near the ferromagnetic phase boundary. 

Jahn-Teller phonon coupling (in Fig.[7J) shows the sim- 
ilar behavior of the local pairing as the HR pairing. 
In contrast to the HR coupling case, Tc is non-zero at 
U = with the pairing mediated by the external cou- 
pling medium of phonons. The results for 2A/W = 0.1 
shows that the phonon pairing-medium provides the at- 
tractive pairing for the small U and as U is increased 
further the superconductivity makes a cross-over to the 
local pairing regime. At large electron-phonon coupling 
{2A/W = 0.15) the metallic region shrunk and showed 
no signs of the cross-over. Given the same coupling 
strength (comparing K and |i^ph| — 2A), the HR system 
has weaker superconductivity than the el-ph systems. 

One significant difference between this model of doubly 
degenerate orbitals (commonly called a,s exE model) and 
the Geo model {t x H model^^,) is that the e x E model 
produces stronger enhancement of Tc in the local pairing 
regime. Authors in Refii^ found no signs for the up-turn 
of Tc near the Mott transition for any studied values of A 
in the Ggo-model, although the enhancement of Tc over 
the conventional theories was evident. It is probably due 
to the easier low spin formation for even fillings, since at 
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odd number of fillings (as 3 in the Cgo model) unpaired 
spins tend to fluctuate via the intersite hopping and make 
the low spin configuration less effective. 

Since the HR and JT couplings induce superconduc- 
tivities separately but in different symmetries, they are 
expected to complete when both are present. For in- 
stance, by turning up the JT coupling constant at a 
fixed finite HR coupling, one can suppress the spin-triplet 
pairing and eventually crosses over to a spin-singlet pair- 
ing, ft can be also anticipated from the mapping of the 
JT coupling to a negative HR hamiltonian in the anti- 
adiabatic limit. FuUerene systems are known to have 
sizable strengths for both interactions, although the JT 
interaction is more dominant for alkali-doped region. A 
recent density functional calculation stud}i2& has sug- 
gested that the fuUerene in the hole-doped region may 
have stronger Hund's rule interaction, with possible ap- 
plications to cation-doped fuUerenes^. 

Finally, we make an observation of an interesting prop- 
erty of the local interaction Eq. (PJ. If we swap the 
spin and orbital indices in the JT coupling, Eq. (PJ, new 
boson fields (pi, couple to electronic spin fluctuations, 
i.e., the coupling terms become X]m['?i("™T ~ "mi) + 
<^2('^.|„|rfmi -l-h.c.)]. Now (pij plays the role of externa/ lo- 
cal spin-fluctuations which mediate the electron pairing. 
All other terms in the hamiltonian are unchanged. The 
ground state now belongs in the spin-triplet and orbital- 
singlet space. The pair susceptibility for the singlet chan- 
nel maps to that of the triplet-channel Eq. (|23|l and all re- 
sults from JT coupling automatically hold for the triplet 
superconductivity. This model provides a robust exam- 
ple for the local spin-fluctuation superconductivity. 



V. CONCLUSION 

We have demonstrated that the local pairing mecha- 
nism in multi-band systems on general grounds. In par- 
ticular, the existence of Hund's rule induced spin-triplet 
s-wave superconductivity is shown. The self-generated 
local multiplet fluctuations in the strong-f7 limit provide 
the pairing medium even with a purely electronic interac- 
tion. This general idea should work for other extensions 
of the model. Further studies on the superconductor- 
ferromagnet phase diagram would clarify the relevance of 
the current model to known experimental systems. Fluc- 
tuation effects beyond the mean field approximation and 
singlet-triplet competition are left for future research. 
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APPENDIX A: ATOMIC LEVEL SCHEME OF 
HUND'S RULE COUPLING 

For 2 electrons on site, there are 6 possible configura- 
tions 

(Al) 

These states are degenerate for 7?muit — 0. -ffmuit — H^r 
of Eq. ^ can be rewritten in the above basis set as 

/ ~2K \ 
-2K 
-K K 
K -K 
K K 
\ K K J 



(A2) 



the solutions of which become 



E = -2K, \ 1/V2{dl^d. 



2i 



d^ d^ 



\vac) 



\vac) 



(A3) 



2K, l/\/2(4|4i + 4T4i)bac) 



As shown in Fig. |21 the ground states are S = 1 spin- 
triplet and the excited states are 5 = spin-singlet 
with orbitally mixed states. The above ground multiplet 
states are written as |Ti), \To) and \Ti) for the triplet 
states (T) of Sz = 1,0,-1, respectively. Similarly, we 
write the above excited states as \Sa), \Sb) and \Sc), re- 
spectively, for spin-singlets (S). With one electron on 
site, the Hund's rule coupling does not lift any degener- 
acy since it is a multi-electron interaction. The same is 
true with 3 electrons, since one can view this in terms of 
the hole interaction. 



APPENDIX B: PERTURBATION THEORY FOR 
PAIRING CORRELATION ENERGY 

The perturbation is defined as an expansion with re- 
spect to the intersite hopping around the interacting 
Hamiltonian in Eq. ||SJ|, 

Ho = UNiN-l)/2 + H^^it (Bl) 

^ = + (B2) 

ma 

First, let us develop a general perturbation theory up to 
the 4-th order. The Heisenberg equation reads 



(B3) 



where the energy eigenvalue and wavefunction are ex- 
panded in terms of the perturbation V as 

E = Eo+Ei+E2+E:i+Ei, (B4) 

IV-) = IV^0)+|V^l) + |^2) + |^3) + |V^4). (B5) 
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The 0-th order terms give 

The n-st order equation becomes 



n—m/ • 



(B6) 



(B7) 



Projecting the unperturbed bra- vector {-tpol to Eq. HB7|) 
for n = 1, one gets the first order correction to the energy 



El = (^olV'IV'o)- 



(B8) 



Introducing an operator Pq which projects out the un- 
perturbed wavefunction li/'o); one can express the wave- 
function to the 1st order as, 



Eq — Hq 



(B9) 



The higher order contributions are computed by induc- 
tion. The useful expressions up to the 4th order are 



1^-2) 



1 



(BIO) 



Ea — Hq 

E2 = (^o\V\^i) = (^iKEo - HoMi) (Bll) 
Ei = (^2\iEo-HoM2)-E2{^M, (B12) 

where we have used a simplification Ei = which is the 
case with our models. 

The unperturbed energies with the Jahn- Teller inter- 
action, Eq. ©, at the total occupation numbers Nd, Nd + 
l,Nd + 2 are 



E{Nd + 2) 


= U- 


8A 


E{Nd + 1) 


= U- 


8A 


E{Nd) 


= U- 


8A 



4tl \M\ 



U AC/2 
7tt 



U 2A[/2 



U 



(B13) 
(B14) 
(B15) 



where A = —8g^/ujp\i is the local Jahn- Teller coupling 
energy. Corresponding wavcfunctions of the leading or- 
der are 

mNd + 2)) = \Sc)^\s,) (B16) 
mNd + l)) = \S,)®cUvac) (B17) 
\tP{Nd)) - \Sc)<E>\vac), (B18) 

with notations given in Appendix IXI 



APPENDIX C: HUBBARD-STRATONOVICH 
DECOUPLING FOR HUND'S RULE COUPLING 

The discrete Hubbard-Stratonovich transforma- 
tioi* (DHST) is given by 

g-Ar(7n„«,3 ^ ^- Ajii +„^) 1 ^ ^Ao-(n„ -r!,3) ^ ^q^^ 



for fermion number operator for a quantum state a 
and cosh A — exp(ATC//2). The continuous Hubbard- 
Stratonovich transformation22(CHST) applies to a gen- 
eral operator A in a gaussian integral 



dx exp(— TTx^ + ^/ttAx). 



(C2) 



For doubly degenerate orbital systems, the interacting 
Hamiltonian Eq. ^ can be rewritten as 



int 



+ K 



E 



(C3) 



E 



E 



'"ma ' '"in' (7 



where N+ = X)™ "™T ^^d N- = J2m'"'ml- Previously, 
each terms in the first two lines have been decoupled by 
the DHST. N^, N^, N+N^ a re ex panded into nma^rn'o-' 
terms before the DHST Eq. ljCl|l is apphed. The terms 
on the third line have been decoupled by 



/ArKalcl^Cy—c^gCs) 



(C4) 



where AT > and no pairs of the indices from {a, /?, 7, 6} 
are the same. When the QMC sampling is incomplete, 
all the symmetry inherent in the original hamiltonian is 
not fully recovered, which is speculated to be a source 
of the sign-problem. This problem can be resolved by a 
different decoupling scheme as follows. We rewrite the 
first line as 



U{Nl , 2 



(C5) 



with Sz = — N^. For the orbital degeneracy Nd — 2, 
the first two terms in the second line in Eq. can be 
written as 

K {ni^nii + n2]n2i - nlTf^2i - n2^nii) 
K 



("-IT - '^2T - "u 



"2i) 

-K^(/V++Af_)-i^(niTn2T+n2imi), (C6) 



K 



where the first term with squared term is decoupled with 
the CHST of a single auxiliary field. This procedure 
preserves the symmetry between orbital 1 and 2. The 
third line in Eq. (|C3|I has been the source of the sign- 
problemSSi. We express this by completing square as 

K{c\^c\^C2iC2^ + 4^4j^CixCi| + • • • ) 



K 

y 

K 



+ y (Af+ -I- iV„) - K{ni^n2i + naxnii). (C7) 
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Similarly to Eq. (|C6|I . the squared term preserves the 
orbital symmetry with a negative coefficient [K > 0). 
With -2(niTn2T +^2X»u ) = -{Nf + N'^) + {N+ + iV_) 
and adding Eqs. I|C5IC7|I . we finally have 

Hint = iK(N+ + N_) + {U- 2K){Nl + N^) 

~ 2"^ - y ('^iT - "2T - nil + ^2i) 
K 

where the A'^i in the first line are expanded into nmaTim'a' 
and decoupled by the DHST, Eq. (EH). Terms in the sec- 
ond and third lines are decoupled by the CHST, Eq. HC2p . 



With the old DHST scheme for all interaction terms, 
there are 6 auxiliary fields for naUp pairs and 4 for 
cJjcjjC^Ci terms, 10 in total per each time slice. The im- 
proved scheme has 2 discrete fields for ni|n2|, Hi|n2| in 
iV^ and 3 continuous fields per time slice. This procedure 
can be extended to higher degeneracy N > 2. Although 
the arrangement in Eq. (|C6p is special for — 2, they 
can be decoupled by the DHST without causing the sign- 
problem. Terms like in Eq. I)C7|I are problematic for the 
sign-problems in the impurity problem. Eq. (jC7|) can be 
readily extended to higher degeneracy and the above pro- 
cedure will likely remove the sign-problem in the Hund's 
rule coupling. 
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